Appendix A 



package analysis; 

5 

import acme.*; 
import java.util.*; 
import java.io.*; 
import java.awt.*; 
10 import java.awt.event.*; 
import javax.swing.*; 

////////////////////////////////////////////////////////////////////////////// 
public class Analysis { 

15 

// Temp for standalone analysis project, rundatastream.java 
public final static short TEMP = 7, OPTICS = 1 * 1024; 

public final static int NORMAL = 0, RAW = 1 , DERIV1 = 2, DERIV2 = 3, 
20 DERIV1 RAW = 4, DERIV2RAW = 5, MELT_OPTICS = 6, 

MELT_TEMPERATURE = 7, MELT_DERIV1 = 8; 
public final static int MAX_CYCLES = 100; 
public final static int MAX_DYES = 4; 
public final static int MAX_SITES = 96; 

25 

// Results 

public final static int PASS = 0; 
public final static int FAIL = 1 ; 

public final static int NO_RESULT = 2; // eg, passive dye 
30 public final static int ND = 3; // Not Determined, IC invalid 



35 



40 



// Dye Types 

public final static int UNUSED = 0; 

public final static int ASSAY = 1 ; 

public final static int I NTERNAL_CONTROL = 2; 

public final static int QIC = 3; 

public final static int PASSIVE = 4; // Historical but needed 

public final static int UNKNOWN = 5; // Qual. Find cone, for this dye 

public final static int STANDARD = 6; // Qual. Dye with known cone. 

// Site Designation 

public final static int SITEJJNKNOWN = 0; 
public final static int SITE_STANDARD = 1 ; 



45 // Data to use 
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public final static int PRIMARY = 0; 

public final static int D2 = 1 ; // 2nd Derivative 

// Analysis Type 
5 public final static int QUALITATIVE = 0; 
public final static int QUANTITATIVE = 1; 

// Threshold mode 

public final static int AUTO THRESH = 0; 
10 public final static int MAN_THRESH = 1 ; 

public static boolean annotate = false; 

// Setup, results... 
M= 15 SitesiteQ; 

Q 

£3 private int analysisType; 

^ // Num Sites 

% 20 private int numSites; 

T II One per dye, site independent 

L. // Primary: 0; 2D: 1 

m int dataTypeQ = new int[MAX_DYES]; 

ifl // Following used for the standards curve, prakash 1/25/00 

Q double dyeSlopeQ = new double[MAX_DYES]; // m: mx+b 

M> double dyeOffsetQ = new double[MAX_DYES]; // b: mx+b 

double linCCQ = new double[MAX_DYES]; 

30 

// standardsLine[0-3][2] 

// Each point is defined by (cycle, Iogb10(concentration)) 

public StdElement standardsLineQD = new StdElement[MAX_DYES][2]; 

public static int stdChannel = 0; 

35 

// IC used: T, IC not used:F 
private boolean useIC; 
private int icDye; 

40 // QIC used: T, QIC not used:F 
private boolean useQIC; 
private int qicDye; 

// Threshold Mode (1 per dye) 
45 private int threshModeQ = new int[MAX_DYES]; 
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// Valid Cycle Number Range for all dyes 

private float validMinCycleQ = new float[MAX_DYES]; 

private float validMaxCycleQ = new f loat[M AX_D YES] ; 

// Cycle Number for noise sub and 3 sigma calculation, 
boolean noise; 

int baselineStartCycleD = new int[MAX_DYES]; 
int baselineEndCycleQ = new int[MAX_DYES]; 

// StdDev baseline for auto threshold detect. User entered, 
// one per dye. 

private double stdDevBaseLineQ = new double[MAX_DYES]; 

// The Max stdDev for a given dye. one per dye 
private float maxStdDevQ = newfloat[MAX_DYES]; 

// This is set to true only if all sites have a valid 

// stdDev. Than only can you calculate the max. 

private boolean maxStdDevValidQ = new boolean[MAX_DYES]; 

// BoxCar Averaging 
private boolean boxCar; 

private int boxCarWidth; // Note Min Value = 2 

// Quantitative Analysis 

public StdElement qtArrQD = new StdElement[MAX_DYES][1]; 



// //////////////////////////////////////////////////////////////////////// 

// Keeps current settings, resets Data (and all calculated values from data) 

II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public void resetData() { 

for(int s = 0; s < numSites; s++) { 
site[s].cycle = 0; 
site[sj. control = false; 
site[s].meltPoints = 0; 

for(int d = 0; d < MAX_DYES; d++) { 
site[s].dye[d].tValid = false; 
site[s].dye[d].tCycle = Of; 
site[s].dye[d].stdDevValid = false; 
site[s].dye[d].slope = 0.; 
site[s].dye[d].offset = 0.; 
site[s].noiseValid[d] = false; 

} 
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} 

// qtArr = null; 

StdElement aQ = new StdElement[1]; 

a[0] = new StdElement(); 

// Site independent 

for(int d = 0; d < MAX_DYES; d++) { 

maxStdDev[d] = Of; 

maxStdDevValid[d] = false; 

qtArr[d] = null; 

qtArr[d] = a; // Reset Quantation 

standardsUne[d][0] = new StdElement(); 
standardsl_ine[d][1] = new StdElement(); 
dyeSlope[d] = 0.; 
dyeOffset[d] = 0.; 
linCC[d] = 0.; 

} 

} 



II IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIHIIIIIIIII 
II Keeps current optics data, redoes all calculations. Eg. may be called 
// after changing Threshold mode from manual to auto. 
II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public void recalcQ { 
int s, cy; 

//System. out.println("Analysis.recalc()"); 

int cQ = new int[numSites]; 

int meltCountQ = new int[numSites]; 

for(s = 0; s < numSites; s++) { 
c[s] = site[s]. cycle; 
meltCount[s] = site[s].meltPoints; 

} 

resetData(); 

for(cy = 0; cy < c[0]; cy++) { 
for(s = 0; s < numSites; s++) { 
addCycle(s, site[s].dye[0].rOptic[cy], site[s].dye[1].rOptic[cy], 
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site[s].dye[2].rOptic[cy], site[s].dye[3].rOptic[cy]); 

} 

} 

} 

5 

ininiiiniiniiiuiiimiuuiumniummuniiiinmuH 

public void setNumSites(int s) { 
if(s <= 0) { 
10 return; 

} 

if(s < numSites) { 
for(int i = s; i < numSites; i++) { 
g 15 site[i] = null; 

1 > 

v"s numSites = s; 

p } 

6 20 

r // niimimiiiuiimiimmiimimniiiiiiiiiimiiii 

¥> public void addCycle(int s, short opO, short op1 , short op2, short op3) { 

ill int c = site[s]. cycle; 

H 25 

, ass. 

//System.out.println("addCycle Site " + s + " OpO " + opO); 

?== site[s].dye[0].rOptic[c] = opO; 

site[s].dye[1].rOptic[c] = op1; 
30 site[s].dye[2].rOptic[c] = op2; 
site[s].dye[3].rOptic[c] = op3; 

site[s].dye[0].pOptic[c] = opO; 
site[s].dye[1].pOptic[c] = op1; 
35 site[s].dye[2].pOptic[c] = op2; 
site[sj.dye[3].pOptic[c] = op3; 

processData(s); 

40 ++site[s].cycle; 
} 

////////////////////////////////////////////////////////////////////////// 

public void addMelt(int s, short sees, int type, short value) { 

45 
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//System.out.println("addMelt Site " + s + " sec " + sees + " type " + type + " 
value " + value); 



- --sz 
■tsss? 



site[s].meltPoints = sees; 

5 

switch(type) { 

//case RunDataStream.OPTICS: 
case OPTICS: 
site[s].mOptic.set(secs, value); 
1 0 site[s] . u pdateM eltDeri v() ; 

break; 

//case RunDataStream.TEMP: 
case TEMP: 

15 site[s].mTemp.set(secs, ((float)value / 1 00f)); 

break; 

} 



BBS 

n 20 



} 



II lllllllllllllllllllllllllllllllllllllllllllltllllllllllllllllllllllillll 
II 0=QI, 1=Qn 

public void setAnalysisType(int a) { 
analysisType = a; 

25 } 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
II To Manually set Threshold limit 
30 // Call this once per dye 

public void setTLimit(int d, float tl) { 
for(int s = 0; s < numSites; s++) { 
site[s].dye[d].tl_imit = tl; 

} 

35 } 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
II For testing quantation only. 
40 // Call this once per dye 

private void setTCycle(int s, int d, float tc) { 

site[s].dye[d].tCycle = tc; 

site[s].dye[d].tValid = true; 

} 

45 
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II llllllllllllllllllllllllllllllllllllllllllllllllilllllllllllllllllllllll 

//0=Auto, 1=Man 

public void setTMode(int d, int tm) { 
threshModefd] = tm; 

5 } 



II llllllllllllllllllllilllllllllllllllllllllllllllllllllllllllllllllllllll 
II Cone, values for Quantitative analysis is set per site per dye 
10 public void setConc(int s, int d, float cone) { 
site[s].dye[d].conc = cone; 

} 



1 5 // llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
//0=Primary, 1=2D 
public void setDataType(int d, int dt) { 
dataType[d] = dt; 

} 

20 



// llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

H // 0=UNKNOWN, 1=STANDARD 

B II In the GUI, SITEJJNKNOWN = 0 and SITE STANDARD = 1 

1* 25 public void setSiteType(int s, int ty) { 

Ci for(int d = 0; d < MAX_DYES; d++) { 

J3 if(!((uselC && d == icDye) || (useQIC && d == qicDye))) { 

'** site[s].dye[d].dyeUsage = ty + 5; 

} 

30 } 
} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
35 // Unused/Std/Passive... 

public void setDyeUsage(int s, int d, int du) { 

switch(du) { 

40 case I NTERNAL_CONTROL: 

for(int si = 0; si < numSites; si++) { 
site[si].dye[d].dyeUsage = du; 

} 

45 useIC = true; 

icDye = d; 
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break; 



case QIC: 

5 for(int si = 0; si < numSites; si++) { 
site[si].dye[d].dyeUsage = du; 

} 

useQIC = true; 
10 qicDye = d; 

break; 

} 

Us ^ 

m 15 

I II lllllllllllllltlltllllllllllllllllllllllllllllllllllllllllllllllllllllll 

Sj // d=Dye, sd = standard dev. Set by User 

=p public void setStdDevbaseline(int d, double sd) { 

C 20 stdDevBaseLine[d] = sd; 



! y // llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

% 25 II IC and Qic 

g public void setlCCycle(int d, int min, int max) { 

jj validMinCyclefd] = (float)min; 

validMaxCycle[d] = (float)max; 

} 

30 

// llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public void setNoiseSubtraction(boolean flag) { 
noise = flag; 

35 } 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public void setBaselineCycle(int dye, int start, int end) { 
40 baselineStartCycle[dye] = start; 
baselineEndCycle[dye] = end; 

} 



45 II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public void setBoxCarAvg(boolean flag, int width) { 
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boxCar = flag; 
boxCarWidth = width; 

} 



// /I//I/I/I/I////////I////////I///J//////I//I///I/I//I//I////I//I///////// 

II Get Thresholds 

public float getTLimit(int s, int d) { 

//System.out.println("Analysis: getTLimit() " + site[s].dye[d].tLimit ); 
10 return site[s].dye[d].tLimit; 

} 



P 

31 

~l 3 

i s 



// //////////////////////////////////////////////////////////////////////// 

15 public float getTCycle(int s, int d) { 

if (site[s].dye[d].tCycle < validMinCycle[d] || site[s].dye[d].tCycle > 
validMaxCycle[d]) 
return Of; 
else 

20 return site[s].dye[d].tCycle; 
} 



II lllllllllllllllllllllllllllllllllllllllllllllllllllllllllltllllillilllll 
25 public float getQICTCycle(int s, int d) { 

int qicDye = getQICDye(); 

float qicTCycle = getTCycle(s, qicDye); 

30 if (useQIC && (qicTCycle > Of)) { 

if (d == qicDye) return qicTCycle; 
return (getTCycle(s,d) / qicTCycle); 

} 

else 

35 return Of; 



} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
40 public boolean getTVaIid(int s, int d) { 
return site[s].dye[d].tValid; 

} 

////////////////////////////////////////////////////////////////////////// 
45 public final double log1 0(double a) { 
if(a>0.){ 
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return (Math.log(a) / Math.log(10.)); 

} 

else { 
return -9.5; 

5 } 
} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
10 public final double Iog10(float a) { 
if(a>0.){ 

return (Math.log((double) a) / Math.log(10.)); 

} 

else { 
15 return -9.5; 

} 

} 



20 // IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIHI 
II Get Results 

II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public int getQLResult(int s, int d) { 

25 int du = site[s].dye[d].dyeUsage; 

// Update IC 

if(uselC &&!site[s].control) { 
updatelC(s); 

30 

} 

if(du == UNUSED || du == PASSIVE) { 
site[s].dye[d].qlResult = NO_RESULT; 

35 

} 

else if(uselC) { 
if(site[s].control) { 

site[s].dye[d].qlResult = site[s].dye[d].tValid ? PASS : FAIL; 

40 } 

else { 

site[s].dye[d].qlResult = ND; 

} 

} 

45 else { 

site[s].dye[d].qlResult = site[s].dye[d].tValid ? PASS : FAIL; 
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} 

return site[s].dye[d].qlResult; 

} 

5 

////////////////////////////////////////////////////////////////////////// 
// Update Internal Control Status 
void updatelC(int s) { 

10 

if(site[s].dye[icDye].tValid) { 

// Also make sure it happened in the specified range 
if((site[s].dye[icDye].tCycle >= validMinCycle[icDye]) && 
15 (site[s].dye[icDye].tCycle <= validMaxCycle[icDye])) { 

site[s]. control = true; 

} 

else { 

20 // Although .tValid, not in the range 
site[s]. control = false; 

} 

} 

else { 

25 site[s]. control = false; 
} 

} 



30 II llllllllllllllllllllllllllllllllilllllllllllllllllllllllllllllllllllllll 
II Update Linear Correlation Coefficient 
// //////////////////////////////////////////////////////////////////////// 
void updateCC(int d) { 

35 double yt, xt; 

double syy = 0., sxy = 0., sxx = 0., ay = 0., ax = 0.; 

if(qtArr[d].length < 2) { 
linCC[d] = 0.; 

40 

return; 

} 

for(int j = 0; j < qtArr[d].length; j++) { 
45 ax += qtArr[d][j].conc; 

ay += qtArr[d][j].avgTCycle; 
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} 

ax /= qtArr[d]. length; 
ay/=qtArr[d].length; 

5 

for(int j = 0; j < qtArr[d].length; j++) { 
xt = qtArr[d]0].conc - ax; 
yt = qtArr[d][j].avgTCycle - ay; 
sxx += xt * xt; 
10 syy += yt * yt; 
sxy += xt * yt; 

} 

linCC[d] = sxy / (Math.sqrt(sxx * syy)); 
M- 15 linCC[d] *= linCC[d]; 



j inmiiuiiMimiiiiiiMMimmimmmiim 

S 20 110. Check for unknown & thresh. 
5 //1. Check IC 

i // 2. Check QIC 

U II 3. Check for at least 2 data points in this qtArr 

nj // 4. Check for unknown to be within knowns 

M= 25 // 5. Sort qtArr and Return unknown cone. Move to addstandard... 

§ II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

p public double getQTResult(int s, int d) { 

double m = 1.0; 

30 

// 0. Check for unknown thresh. 

if(!site[s].dye[d].tValid || (site[s].dye[d].dyeUsage != UNKNOWN)) { 
return 0.; 

} 

35 

II A. Check IC 
if(uselC) { 
if(!site[s].dye[icDye].tValid) { 
return 0.; 

40 } 
} 

// 2. Check QIC 
//todo prakash. 

45 // Should wait for all thresholds/site before constructing qtArr. 
if(useQIC) { 



AttyRef#22660-0025B US 



if(!site[s].dye[qicDye].tValid){ 
return 0.; 

} 

else { 

5 m = 1 . / site[s].dye[qicDye].tCycle; 
} 

} 

// 3. Check for at least 2 data points in this qtArr 
10 if(qtArr[d].length < 2) { 
return 0.; 

} 

m site[s].dye[d].conc = (float) Math.pow(10., (dyeSlope[d] * 

O 15 (site[s].dye[d].tCycle * m) + dyeOffset[d])); 

fU // 4. Check for the cone to be within .5 Log 

if( (Iog10(site[s].dye[d].conc) > standardsl_ine[d][0].conc) || 
J (Iog10(site[s].dye[d].conc) < standardsLine[d][1].conc)) { 

% 20 site[s].dye[d].conc = Of; 

?. return site[s].dye[d].conc; 

fy > 

S 25 

Q II lllllllllllllllllllllllllllllllllllllllllllllllllllllllllillllllllilllll 

%a II Sort the elements in the Quantation Array, 

void sort(StdElement aQ) { 

30 boolean done; 

StdElement se = new StdElement(); 

if(a. length < 2) { 
return; 

35 } 

do{ 
done = true; 

40 for(int j = 0; j < (a.length - 1 ); j++) { 

if(a[j].avgTCycle > a[j + 1].avgTCycle) { 

done = false; 

se = aG]; 

aD] = aO + 1]; 
45 a0 + 1] = se; 
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break; 

} 

} 

} 

5 while(ldone); 
} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
10 // Sort the elements in the Melt Peaks Array, 
void sort(MeltElement meltElementsArrayO) { 

boolean done; 

MeltElement me = new MeltElement(); 

g //Debug.log ("sort: MeltElement array with " + meltElementsArray. length); 

Jri if(meltElementsArray.length < 2) { 

return; 

n i 

•6 20 

5 do { 

done = true; 

HI for(int j = 0; j < (meltElementsArray.length - 1 ); j++) { 

M 25 if(meltElementsArray[j].d1 Peak > meltElementsArrayO + 1].d1 Peak) { 
jg done = false; 

O me = meltElementsArrayO]; 

^ meltElementsArrayO] = meltElementsArrayO + 1]; 

meltElementsArrayO + 1] = me; 

30 

break; 

} 

} 

} 

35 while(ldone); 
} 



iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiniiiiiiiiiiim 

40 // Update data used for drawing the Line fit to standards. 

// 

// standardsLine is similar to qtArrQ but adds 2 points, one at 
// cone +.5(log) and the other at cone -.5 (log). 

II llllllllllllllllllllllllllllllllllllllllllllllllllllllllltllllllllllllll 
45 void updateStandards(int d) { 
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int e = qtArr[d]. length - 1 ; 

double cone = qtArr[d][e].conc - .5; 

standardsLine[d][0].conc = qtArr[d][0].conc + .5; 
5 standardsLine[d][0].avgTCycle = (standardsLine[d][0].conc - dyeOffset[d]) 

/ dyeSlope[d]; 

if(conc > 0.) { 
standardsLine[d][1].conc = cone; 
10 standardsLine[d][1].avgTCycle = (cone - dyeOffsetfd]) / dyeSlope[d]; 
} 

else { 

standardsLine[d][1].conc = 0.; 

standardsLine[d][1].avgTCycle = (-1 * dyeOffset[d] / dyeSlope[d]); 

15 } 
} 



N II lllllllllllllllllilllllllllllllllllllllllllllillilllllllllllllllllllllll 

± 20 // Get Control Result (Pass/Fail) 

5 II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

public boolean getControl(int s, int d) { 
L return site[s].control; 

m } 



ii iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii 

public float getConc(int s, int d) { 
return site[s].dye[d].conc; 

30 } 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public int getDyeUsage(int s, int d) { 
35 return site[s].dye[d].dyeUsage; 
} 



// llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllil 
40 public double getDyeSlope() { 
return dyeSlope[stdChannel]; 

} 



45 II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public double getDyeOffset() { 
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return dyeOffset[stdChannel]; 

} 



5 II lllllllllllllllllllllllllllllllllllllllllllllllllllllllllillllllllllllll 
II Linear Correlation Coefficient 
public double getCC() { 
updateCC(stdChannel); 

10 return linCC[stdChannel]; 
} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
15 public float getAnaData(int dataType, int s, int d, int c) { 

float retVal = Of; 

if (c < 0) c=0; 

switch (dataType) { 

case NORMAL: 

if (c >=site[s]. cycle) c=site[s]. cycle - 1; 
if(d < 4 && d >= 0) { 
retVal = site[s].dye[d].pOptic[c]; 

} 

break; 

30 case DERIV1 : 
break; 

case DERIV2: 
if (c >=site[s]. cycle) c=site[s].cycle - 1 ; 
35 if(d < 4 && d >= 0) { 

retVal = site[s].dye[d].d2pOptic[c]; 

} 

break; 

40 case MELT_DERIV1: 

if (c >=site[sJ.meltPoints) c=site[s].meltPoints - 1; 
if(c < site[s].meltPoints && c>= 0) { 
retVal = site[s].d1mOptic.get(c); 

} 

45 break; 




% 20 
2 25 
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case MELT_OPTICS: 
if (c >=site[s].meltPoints) c=site[s].meltPoints - 1; 
if(c < site[s].meltPoints && c >= 0) { 
retVal = site[s].mOptic.get(c); 

5 } 

break; 

case MELT_TEMPERATURE: 
if (c >=site[s].meltPoints) c=site[s].meltPoints - 1 ; 
10 if(c < site[s].meltPoints && c >= 0) { 

retVal = site[s].mTemp.get(c); 

} 

break; 

} 

S return retVal; 

y } 

l 78? 

J 20 ///////////////////////////////» 
45 public int getlCDye() { 

n return icDye; 

ft } 

h 25 // uimimiiiimuuimumimiuuHmmiiiuimmmium 

jf public boolean iCEnabled() { 

:-• return useIC; 



30 

II llllllllilllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
II Returns the temp assoc. with the Melt Peak, 
public double getMeltTemp(int s, int index) { 
return site[s].getMeltTemp(index); 

35 } 



II llllllllllllltllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
II Returns the Melt Limit. Peak value reported only when greater. 
40 public double getMeltl_imit(int s) { 
return site[s].meltPeakl_imit; 

} 



45 // imuiiimmiiniiiiiiiiiiuininuiinmiiuiiimuiiuii 

II Returns the temp assoc. with the Melt Peak. 
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public int getMeltCount(int s) { 
if (s>0 && s<numSites) 

return site[s].getMeltPeakCount(); 
else 

5 return 0; 
} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
10 public int getQICDye() { 
return qicDye; 

} 



£ 15 // llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllilllllllll 
% public boolean qicEnabled() { 

m, return useQIC; 

,p II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

« public int getTMode(int d) { 

h* return threshMode[d]; 

W } 
5 25 

P // lllllllllllllllllllllllllllllllllllllllllllllllllllllllllll^ 

int getlCStartCycleO { 
return (int)validMinCycle[icDye]; 

30 } 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
int getlCEndCycle() { 
35 return (int)validMaxCycle[icDye]; 
} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
40 void processData(int s) { 

if(boxCar) { 
boxCarAvg(s); 

} 

45 

if(noise) { 
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removeNoise(s); 

} 

updateThresholds(s); 

// Update qtArr's. Do quantation when results are requested. 
if(analysisType == QUANTITATIVE) 
u pdateQua ntitative(s); 

} 



II Apply this to raw Data 
void boxCarAvg(int s) { 

float sum; 
int i; 

if(site[s].cycle < 1){ 
return; 

} 

if(site[s].cycle + 1 >= boxCarWidth && boxCarWidth > 1) { 

for(int d = 0; d < MAX_DYES; d++) { 
sum = Of; 

for(i = (site[s].cycle + 1 - boxCarWidth); i < site[s].cycle + 1; i++) { 
sum += site[s].dye[d].rOptic[i]; 

} 

site[s].dye[d].pOptic[site[s].cycle] = sum / boxCarWidth; 

} 

} 

} 



llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
void removeNoise(int s) { 

int c = site[s].cycle; 
float temp; 



for(int d = 0; d < MAX_DYES; d++) { 
if(c >= (baselineEndCycle[d] - 1)) { 
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if(site[s].noiseValid[d]) { 
site[s].dye[d].pOptic[c] -= (site[s].dye[d].slope * c + site[s].dye[d].offset); 
site[s].dye[d].pOptic[c] -= site[s].dye[d].noiseAvg; 

//if (s==0 && d==0) { 

// Logger.log("Cycle "+c+ " slope "+site[s].dye[d]. slope + 

// " offset " + site[s].dye[d]. offset + " pOptic " + site[s].dye[d].pOptic[c]); 

//} 

} 

else { 
temp = Of; 

// Calculate Average noise 

baselineStartCycle[d] = (baselineStartCycle[d] < 1 ) ? 1 : 
baselineStartCycle[d]; 

site[s].dye[d]. slope = 0.; 
site[s].dye[d].offset = 0.; 

site[s].dye[d].leastSquaresLineFit(baselineStartCycle[d]-1, 
baselineEndCycle[d]-1 ); 

for(int i = 0; i <= (baselineEndCycle[d] - 1 ); { 

site[s].dye[d].pOptic[i] -= (site[s].dye[d].slope * i + site[s].dye[d]. offset); 



for(int i=baselineStartCycle[d]-1 ; i<=baselineEndCycle[d]-1 ; i++) { 
temp = temp + site[s].dye[d].pOptic[i]; 

} 

site[s].dye[d].noiseAvg = temp / (baselineEndCycle[d] - 
baselineStartCycle[d] + 1); 

// Remove noise 

for(int i=0; i <= (baselineEndCycle[d]-1 ); i++) { 
site[s].dye[d].pOptic[i] -= site[s].dye[d].noiseAvg; 

} 

site[s].noiseValid[d] = true; 

} 

} 

} 

} 



II llllllilllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
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void updateThresholds(int s) { 

for(int d = 0; d < MAX_DYES; d++) { 

5 // Update Derivative 
update2D(s, d); 

if(dataType[d] == PRIMARY) { 
if(threshMode[d] == MAN_THRESH) { 
10 updateThreshPDMan(s, d); 

} 

else { 

updateThreshPDAuto(s, d); 

} 

15 } 

else{ 

if(threshMode[d] == MAN_THRESH) { 
updateThresh2DMan(s, d); 

} 

20 else { 

updateThresh2DAuto(s, d); 

} 

} 

} 

25 } 



II llllllllllllllllllllllllilllllllllllllllllllllllllllllllllllllllllllllll 
int updateThreshPDMan(int s, int d) { 

30 

int c = site[s]. cycle; 

int du = site[s].dye[d].dyeUsage; 

if(du == UNUSED || du == PASSIVE) { 
35 return 0; 

} 

if(noise) { 

if(c <= baselineEndCycle[d]) { 
40 return 0; 
} 

} 

if(!site[s].dye[d].tValid) { 
45 if(site[s].dye[d].pOptic[c] >= site[s].dye[d].tl_imit) { 
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// Optic exceeded limit, calculate cycle 
if(c>=1){ 

site[s].dye[d].tValid = true; 

5 LinearFit I; 

I = new LinearFit(c - 1, site[s].dye[d].pOptic[c - 1], c, 
site[s].dye[d].pOptic[c]); 

10 //zero based 

site[s].dye[d].tCycle = l.fitY(site[s].dye[d].tLimit) + 1f; 

} 

} 

} 

15 return 0; 
} 



II llllllllllllillllllllillllllllllllllllllllllllllllllllllllllllllllllllll 
20 // When not to find the Threshold crossing: 
// 

// 1 . Unused Dye 
// 2. Passive dye 
// 3. Already found (.tValid) 
25 // 4. Not enough cycles (2D) 

// 5. All dyes don't have valid stdDev Auto 
IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIH 
int updateThreshPDAuto(int s, int d) { 

30 int c = site[s].cycle; 
float sum, temp; 
int du = site[s].dye[d].dyeUsage; 

if(du == UNUSED || du == PASSIVE) { 
35 return 0; 

} 

if(c <= baselineEndCycle[d]) { 
return 0; 

40 } 

if(maxStdDevValid[d] &&!site[s].dye[d].tValid) { 

// Look for signal crossing 
45 if(site[s].dye[d].pOptic[c] > site[s].dye[d].tl_imit) { 



Atty Ref# 22660-0025B US 



144 



LinearFit I; 

I = new LinearFit(c - 1, site[s].dye[d].pOptic[c - 1], c, site[s].dye[d].pOptic[c]); 
// Add one to match graph 

site[s].dye[d].tCycle = LfitY(site[s].dye[d].tLimit) + 1 .Of; 
site[s].dye[d].tValid = true; 

} 

} 

else if(!maxStdDevValid[d] &&!site[s].dye[d].tValid) { 

// If enough data, calculate stdDev 
// No need to check crossing yet. 
if(c >= baselineEndCycle[d]) { 
if((baselineEndCycle[d] - baselineStartCycle[d]) > 1) { 

// mean 
sum = Of; 

for(c = (baselineStartCycle[d] - 1); c <= (baselineEndCycle[d] - 1); C++) { 
sum = sum + site[s].dye[d].pOptic[c]; 

} 

site[s].dye[d].mean = sum / (baselineEndCycle[d] - baselineStartCycle[d] + 



// stdDev 
sum = Of; 

for(c = (baselineStartCyclefd] - 1); c <= (baselineEndCycle[d] - 1); C++) { 
temp = site[s].dye[d].pOptic[c] - site[s].dye[d].mean; 
sum = sum + temp * temp; 

} 

site[s].dye[d].stdDev = (float) Math.sqrtfsum / (baselineEndCyclefd] - 
baselineStartCycle[d])); 

site[s].dye[d].stdDevValid = true; 

setMaxStdDev(d); 

} 

} 

} 

return 0; 

} 
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// llillllllllllllllllllillllllllillllllllilllllllllllillllilllllllll^ 
II This function calculates the Cycle Threshold for Primary Data with 
// a manual threshold limit set by the user. 
5 // IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII^ 
int updateThresh2DMan(int s, int d) { 

int du = site[s].dye[d].dyeUsage; 

10 // Because the calculation for D2 is lagging 2 cycles back, 
int c = site[s]. cycle - 2; 

if(du == UNUSED || du == PASSIVE) { 
return 0; 

£ 15 } 

-SMS:* 

fU if(c < 6) { 

-\t return 0; 

i } 

Q 20 

*F if(noise) { 

f if(c <= baselineEndCycle[d]) { 

return 0; 

FU } 

fS 25 } 

JJ // Look for peak 

// When c == 6, Possible valid D2's are at c2(c-4), c3(c-3), c4(c-2) 
if((site[s].dye[d].d2pOptic[c - 3] > site[s].dye[d].d2pOptic[c - 4]) && 
30 (site[s].dye[d].d2pOptic[c - 3] >= site[s].dye[d].d2pOptic[c - 2])) { 

PeakFinder peakFinder = new PeakFinder((float) (c - 4), 
site[s].dye[d].d2pOptic[c - 4], 

(float) (c - 3), site[s].dye[d].d2pOptic[c - 3], (float) (c - 2), 
35 site[s].dye[d].d2pOptic[c - 2]); 

// Look for signal crossing 

if (peakFinder, peak > site[s].dye[d].tLimit) { 

40 // peak exceeded limit, calculate cycle 

// Note: peak is 3 cycles back from here 
if(site[s].dye[d].tValid) { 

if (site[s].dye[d].tCycle < peakFinder.cycle + 1 .Of) { 
45 site[s].dye[d].tCycle = peakFinder.cycle + 1 .Of; 

} 
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} 

else { 

site[s].dye[d].tValid = true; 
5 site[s].dye[d].tCycle = peakFinder.cycle + 1 .Of; 

} 

} 

} 

return 0; 

10 } 



ii iiiimiimniiiiiimiimnmiiiiiiimiiiiiiiiiuiiiiiiiniii 

J; 15 int updateThresh2DAuto(int s, int d) { 

^1 int du = site[s].dye[d].dyeUsage; 

\|| float sum, temp; 

£ 'nt cy; 

0 20 

£ II Because the calculation for D2 is lagging 2 cycles back. 

■ int c = site[s]. cycle - 2; 

W if(du == UNUSED || du == PASSIVE) { 

^ 25 return 0; 



if(c < 6) { 
return 0; 

30 } 

if(c <= baselineEndCycle[d]) { 
return 0; 

} 

35 

if(maxStdDevValid[d]) { 

// Look for signal crossing, ie Look for peak 
// When c == 6, Possible valid D2's are at c2(c-4), c3(c-3), c4(c-2) 
40 if(c < (baselineEndCycle[d] + 3)) { 
return 0; 

} 

if((site[s].dye[d].d2pOptic[c - 3] >= site[s].dye[d].d2pOptic[c - 4]) && 
45 (site[s].dye[d].d2pOptic[c - 3] > site[s].dye[d].d2pOptic[c - 2])) { 
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PeakFinder m = new PeakFinder((float) (c - 4), site[s].dye[d].d2pOptic[c - 4], 
(float) (c - 3), site[s].dye[d].d2p0ptic[c - 3], (float) (c - 2), 
site[s].dye[d].d2pOptic[c - 2]); 

5 // Look for signal crossing 

if(m.peak > site[s].dye[d].tl_imit) { 

if (site[s].dye[d].tValid) { 
if (site[s].dye[d].tCycle < m.cycle + 1f) { 
10 site[s].dye[d].tCycle = m.cycle + 1f; 

} 

} 

else{ 

// peak exceeded limit, calculate cycle 
15 site[s].dye[d].tValid = true; 

site[s].dye[d].tCycle = m.cycle + 1f; 

} 

} 

} 

20 } 

else if(!maxStdDevValid[d] &&!site[s].dye[d].tValid) { 

// If enough data, calculate stdDev 
// No need to check crossing yet. 
25 if(c >= baselineEndCycle[d]) { 

if((baselineEndCycle[d] - baselineStartCycle[d]) > 1) { 

// mean 
sum = Of; 

30 

for(c = (baselineStartCycle[d] - 1); c <= (baselineEndCycle[d] - 1); C++) { 
sum = sum + site[s].dye[d].d2pOptic[c]; 

} 

35 // Changed 1/12/00 as per SCR 129. 

// sum = sum + site[s].dye[d].pOptic[c]; 

site[s].dye[d].mean = sum / (baselineEndCycle[dl - baselineStartCycle[dl + 

1); 

40 

// stdDev 
sum = Of; 

for(c = (baselineStartCycle[d] - 1); c <= (baselineEndCyclefdl - 1); C++) { 

45 

// Changed 1/12/00 as per SCR 129. 
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// temp = site[s].dye[d].pOptic[c] - site[s].dye[d].mean; 

temp = site[s].dye[d].d2pOptic[c] - site[s].dye[d].mean; 
sum = sum + temp * temp; 

5 } 

site[s].dye[d].stdDev = (float) Math.sqrt(sum / (baselineEndCycle[d] - 
baselineStartCycle[d])); 

site[s].dye[d].stdDevValid = true; 

10 

setMaxStdDev(d); 

} 

} 

} 

U 15 

O return 0; 

9 > 



.J 20 // IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIHIIIIIIIIIIIIIHin 

H II Update 2nd Deriv for optic data 

T II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

^ void update2D(int s, int d) { 

11 25 int c = site[s].cycle; 
ks float mult = 6.25f; 

N= if (c<4) 

return; 

30 

//D2 

if(c < MAX_CYCLES - 1 && c> 2) { 
// n=3 thru n-2 

35 

/* 

//float mult = 5f; 

site[s].dye[d].d2pOptic[c - 2] = (site[s].dye[d].arD1Dye[c - 1] - 

site[s].dye[d].arD1 Dye[c - 3]) / 2f * mult; 
40 site[s].dye[d].d2pOptic[c - 1] = (site[s].dye[d].arD1 Dye[c] - 

site[s].dye[d].arD1 Dye[c - 2]) / 2f * mult; 
site[s].dye[d].d2pOptic[c] = (site[s].dye[d].arD1Dye[c] - 
site[s].dye[d].arD1Dye[c - 1]) * mult; 

7 

45 site[s].dye[d].d2pOptic[c-2] = (site[s].dye[d].pOptic[c] - 

2f * site[s].dye[d].pOptic[c-2] + 
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1=5 



site[s].dye[d].pOptic[c-4]) * mult; 

site[s].dye[d].d2pOptic[c-1] = (2f * site[s].dye[d].pOptic[c] 
3f * site[s].dye[d].pOptic[c-1] + 
5 site[s].dye[d].pOptic[c-3]) * mult; 

site[s].dye[d].d2pOptic[c] = (site[s].dye[d].pOptic[c] - 
2f * site[s].dye[d].pOptic[c-1] + 
site[s].dye[d].pOptic[c-2]) * 2 * mult; 

10 } 

else{ 

site[s].dye[d].d2pOptic[c] = Of; 

} 



M= 15 

U 



} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
II Update qtArfs (1 per dye - site independent). 
// Only if std: only with valid thresh 
20 II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
void updateQuantitative(int s) { 

for(int d = 0; d < MAX_DYES; d++) { 
if(site[s].dye[d].dyellsage == STANDARD) { 
25 // if(site[s].dye[d].tValid) { 

if( (useQIC && (getTCycle(s, qicDye) > Of)) || getTCycle(s, d) > Of ) { 
addStandard(s, d); 
//updateStandards(d); 
30 LeastSquares Is = new LeastSquares(qtArr[d], d); 

dyeSlope[d] = ls.getSlope(); 
dyeOffsetfd] = ls.getOffset(); 
updateStandards(d); 

} 

35 } 



} 



} 



40 II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
II 

II Add a std Element to the qIArr if appropriate. 
//If QIC used -valid 
// If IC used - valid 
45 // Sort if more than 1 element 

II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
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int addStandard(int s, int d) { 
int i; 

float tCycle; 

if(!site[s].dye[d].tValid || getTCycle(s,d) <= Of ) { 
return 0; 

} 

if(site[s].dye[d].conc < 1 0E-5f) { 
return 0; 

} 

if (useQIC) { 
tCycle = getQICTCycle(s,d); 

} 

else { 

tCycle = getTCycle(s,d); 

} 

if(qtArr[d][0].conc<-9){ 
// Initialise 

qtArr[d][0].conc = log 10(site[s].dye[d]. cone); 
qtArr[d][0].avgTCycle = tCycle; 
qtArr[d][0].nElements = 1; 
return 0; 

} 

else { 

// Look for cone in array 

for(i = 0; i < qtArrfd]. length; i++) { 

if(Math.abs(qtArr[d][i].conc - Iog10(site[s].dye[d].conc)) < .05) { 
qtArr[d][i].avgTCycle = ((qtArr[d][i].avgTCycle * qtArr[d][i].nElements) + 

tCycle) / (qtArr[d][i].nElements + 1 ); 
qtArr[d][i].nElements += 1; 

// May need to be resorted 
if(qtArr[d].length > 1){ 
sort(qtArr[d]); 

} 

return 0; 

} 

} 

// Cone not found, add new element to array 
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StdElement tempArrQ = new StdElement[qtArr[d].length + 1]; 

// Initialise tempArr 
for(i = 0; i < tempArr.length; i++) { 
5 tempArr[i] = new StdElement(); 
} 

System.arraycopy(qtArr[d], 0, tempArr, 0, qtArr[d].length); 

10 tempArrftempArr.length - 1].conc = Iog10(site[s].dye[d].conc); 
tempArr[tempArr.length - 1].avgTCycle = tCycle; 
tempArr[tempArr.length - 1].nElements = 1; 
qtArr[d] = tempArr; 

15 //Sort 

sort(qtArr[d]); 

} 

return 0; 

20 } 



II IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIHII 
void setMaxStdDev(int d) { 

25 

maxStdDevValid[d] = true; 
int s; 

30 maxStdDev[d] = Of; 

for(s = 0; s < numSites; s++) { 
if(site[s].dye[d].stdDevValid) { 
if(site[s].dye[d].stdDev > maxStdDev[d]) { 
35 maxStdDev[d] = site[s].dye[d].stdDev; 

} 

} 

else { 

maxStdDevValid[d] = false; 
40 maxStdDev[d] = Of; 

return; 

} 

} 

45 

if(maxStdDevValid[d]) { 
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// All sites have stdDevValid for dye d, 
// Calculate Threshold limits 
for(s = 0; s < numSites; s++) { 

site[s].dye[d].tLimit = (float)(stdDevBasel_ine[d] * maxStdDev[d]); 

//System.out.println("stdDevBaseLine[d] " + stdDevBasel_ine[d] + 

// "maxStdDev[d] " + maxStdDev[d] + 

// " setMaxStdDev " + site[s].dye[d].tLimit ); 

} 

} 



II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 
public Analysis() { 
this(MAX_SITES); 

} 

public Analysis(int ns) { 

numSites = ns; 

site = new Site[numSites]; 

for(int i = 0; i < numSites; i++) { 
site[i] = new Analysis.Site(); 

} 

analysisType = QUALITATIVE; 

useQIC = false; 
qicDye = 0; 
useIC = false; 
icDye = 0; 

boxCar = false; 
boxCarWidth = 0; 

// Default to match noise sub with primary data. 
// noise = false; 

for(int i = 0; i < MAX_DYES; i++) { 
threshMode[i] = AUTO_THRESH; 
stdDevBasel_ine[i] = 5f; 
maxStdDev[i] = Of; 
maxStdDevValidp] = false; 
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dataType[i] = PRIMARY; 
qtArr[i][0] = new StdElement(); 
baselineStartCycle[i] = 3; 
baselineEndCycle[i] = 8; 

5 

//Standards Curve, prakash 1/25/00 
standardsLine[i][0] = new StdElementQ; 
standardsLine[i][1] = new StdElement(); 

10 // Optics must cross threshold in this range 
validMinCycle[i] = 3f; 
validMaxCycle[i] = 60f; 

} 

} 



IJ class Site { 

5 20 Dye dyeQ = new Dye[MAX_DYES]; 
= // Melt Peak Analysis 

§=* private Array.Short mOptic = new Array.Short(32); 

TU private Array.Float mTemp = new Array.Float(32); 

H" 25 private Array.Float d1 mOptic = new Array. Float(32); 
^2 private MeltElement mPeaksQ = new MeltElement[1]; 

// Possible to set per site in future, 
private double meltPeakLimit = 10.; 

30 

// Melt peaks processed 
private boolean meltPeaksValid; 

// Current Cycle Number 
35 int cycle; 

// Number of MeltData points 
private int meltPoints; 

40 // IC/QIC passed:T; failediF 
boolean control; 

// Noise 

boolean noiseValidQ = new boolean [MAX_DYES]; 

45 

SiteQ { 
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// Initialise dyes 

for(int i = 0; i < MAX_DYES; i++) { 
dye[i] = new Dye(); 
noiseValid[i] = false; 

} 

cycle = 0; 
meltPoints = 0; 
meltPeaksValid = false; 
control = false; 

mPeaks[0] = new MeltElement(); 

} 

private void updateMeltDeriv() { 
meltPeaksValid = false; 

if (meltPoints < 1){ 
d1mOptic.set(0, Of); 

} 

else if(meltPoints == 1 ) { 

d1mOptic.set(1, (mOptic.get(1 ) - mOptic.get(O)) * -5f); 

} 

else { 

// Recalc the 2nd last value, and the last value 
d1mOptic.set(meltPoints-1, (mOptic.get(meltPoints) - 

mOptic.get(meltPoints-2)) / 2f * -5f); 

d1mOptic.set(meltPoints, (mOptic.get(meltPoints) - 

mOptic.get(meltPoints-l)) * -5f); 
} 

} 

// Return number of Melt Peaks detected, 
private int getMeltPeakCount() { 

if (ImeltPeaksValid) 
detectMeltPeaks(); 

return (mPeaks[0].temp < 0.) ? 0 : mPeaks.length; 

} 

// Return number of Melt Temp Associated with Peak, 
private double getMeltTemp(int index) { 
if (index < getMeltPeakCount()) 
return mPeaks[index].temp; 
else 
return Of; 
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} 

// Find all peaks in 1st Deriv of Melt Optic 
private void detectMeltPeaks() { 

5 

if (meltPoints < 2) return; 

if (ImeltPeaksValid) { 
meltPeaksValid = true; 
10 mPeaks = new MeltElement[1]; 

mPeaks[0] = new MeltElement(); 

// Debug.log("detectMP, length " + mPeaks.length); 



15 



for (int i=1 ; i<meltPoints-1; i++) { 

if( ( dlmOptic.get(i) > d1mOptic.get(i-1) ) && 
( dlmOptic.get(i) >= d1mOptic.get(i+1) ) ) { 



PeakFinder peakFinder = new PeakFinder((float)(i-1), 
20 (float)d1mOptic.get(i-1), 

(float)(i), (float)dlmOptic.get(i), (float)(i+1), 
(float)d1mOptic.get(i+1)); 

// Look for signal crossing 
25 if(peakFinder.peak > meltPeakLimit) { 



30 in this case. 



if (mPeaks[0].temp < 0.) { 

mPeaks[0].d1Peak = peakFinder. peak; 
mPeaks[0].temp = mTemp.get(O) + peakFinder.cycle; // Temp, 

} 

else { 

MeltElement tempAQ = new MeltElement[mPeaks.length+1]; 

35 // Initialise tempA 

for(int j = 0; j < tempA.length; j++) { 
tempA[j] = new MeltElement(); 

} 

40 System.arraycopy(mPeaks, 0, tempA, 0, mPeaks.length); 

tempA[tempA.length-1].d1Peak = peakFinder.peak; 
tempA[tempA.length-1].temp = mTemp.get(O) + 
peakFinder.cycle; // Temp, in this case. 
45 mPeaks = tempA; 

} 
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} 

} 

} 

} 

//Debug. log(" detectMeltPeaks() mPeaks. length " + mPeaks. length); 
if (mPeaks. length > 1) 
sort(mPeaks); 

} 

} 

// iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiim 

class Dye { 
// Data Arrays 

short rOpticQ = new short[MAX_CYCLES]; 
float pOpticQ = newfloat[MAX_CYCLES]; 

// 2nd derivative 

float d2p0pticQ = new float[MAX_CYCLES]; 

// Threshold limit 
float tLimit; 
float tCycle; 

// Indicates if signal crossed the Threshold Limit 
boolean tValid; 

// Qualitative Result 
int qlResult; 

//IC, QIC, Unused,... 
int dyeUsage; 

// true = Std; false = Unkn 
boolean std; 

// Dye Concentration 
float cone; 

// Background Noise Value 
float noiseAvg; 

// Std Dev, Mean calculated, one per dye per site 
boolean stdDevValid; 
float std Dev; 
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float mean; 



// For slope removal. One per dye per site 
double slope; 
5 double offset; 

Dye() { 

// Initialise arrays 
10 for(int i = 0; i < MAX_CYCLES; i++) { 
rOptic[i] = 0; 
pOptic[i] = Of; 
d2pOptic[i] = Of; 

} 

g 15 

X // Default Man Threshold, dyeUsage, tValid 

5 qlResult = 0; 

sj tLimit = 200f; 

i tCycle = Of; 

Q 20 tValid = false; 

=p dyeUsage = ASSAY; 

■* std = false; 

Jf; conc = 10E-6f; 

j y noiseAvg = Of; 

£ 25 stdDevValid = false; 

£j stdDev = Of; 

ft mean = Of; 

slope = 0.; 

offset = 0.; 

30 } 

void endPointl_ineFit(int start, int end) { 

slope = (pOptic[end] - pOptic[start]) / (double)(end - start); 

35 if ((slope * end) != 0.) { 

offset = pOptic[end] / (slope * end); 

} 

else { 
offset = 0.; 

40 } 
} 

void leastSquaresLineFit(int start, int end) { 

45 if ((end - start) < 2) { 

return; 
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} 

LeastSquares Is = new LeastSquares(pOptic, start, end); 
slope = ls.getSlope(); 

5 if ((slope * end) != 0.) { 
offset = ls.getOffset(); 

} 

else { 
offset = 0.; 

10 } 
} 

} 



£15 // llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllilllllllll 

"g public class StdElement { 

n\ public double cone; 

nj public double avgTCycle; 

£ int nElements; 

6 20 

£ StdElement() { 

; cone = -10.; 

H avgTCycle = 0.; 

nElements = 0; 

25 } 

*? } 

// //////////////////////////////////////////////////////////////////////// 

public class MeltElement { 
30 public double temp = -1.; 

public double d1 Peak = -1 .; 

} 



35 

illlllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

public static void main(String argsQ) { 

40 int s, d, c, cy; 

Analysis a = new Analysis(); 

// For reading data from Excel 
Vector vFam = new Vector(16); 
45 vFam.setSize(16); 

Vector vTet = new Vector(1 6); 
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vTet.setSize(16); 

Vector vTam = new Vector(16); 

vTam.setSize(1 6); 

Vector vRox = new Vector(16); 

vRox.setSize(16); 

// Analysis Type 

a.setAnalysisType(QUALITATIVE); 
//a.setAnalysisType(QUANTITATIVE); 

a.setNumSites(16); 



for (d=0; d<MAX_DYES; d++) { 

//a.setDataType(d, D2); // Set Up Data Type 
a.setDataType(d, PRIMARY); 

a.threshMode[d] = AUTO_THRESH; // Set Thresh Mode 
//a.threshMode[d] = MANTHRESH; 

a.stdDevBaseLine[d] = 5.; 
} 

// Set Threshold 
//a.setTLimit(0, 10f); 
//a.setTLimit(1, 10f); 
//a.setTLimit(2, 10f); 
//a.setTLimit(3, 10f); 

// Test BoxCar Avg 
a.setBoxCarAvg(true, 3); 

// Test QIC Dye 
a.setDyeUsage(0, 1, QIC); 

// Test Background Noise Subtraction 
a.setNoiseSubtraction(true); 

// Valid Min, Max Cycle defaults to 3, 60 
//a.setlCCycle(3, 30, 60); 

// Add Data Thresholds and cycle crossings are calculated as soon as 
// enough data has accumulated. 
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try { 

Buffered Reader in = new BufferedReader(new 
FileReader("data5.csv")); 

String str; 

// Throw away first 2 lines 
str = in.readl_ine(); 
str = in.readl_ine(); 

while ((str = in.readl_ine()) != null) { 
//Debug.log(str.length()+" "+ str); 
StringTokenizer t = new StringTokenizer(str, ","); 

for (int i=0; i<16; i++) 
if (t.hasMoreTokens()) 
vFam.setElementAt( (lnteger.valueOf(t.nextToken())) ( i); 

for (int i=0; i<16; i++) 
if (t.hasMoreTokens()) 
vTet.setElementAt((lnteger.valueOf(t.nextToken() )), i ); 

for (int i=0; i<16; i++) 
if (t.hasMoreTokens()) 
vTam.setElementAt((lnteger.valueOf(t.nextToken() )), i ); 

for (int i=0; i<16; i++) 
if (t.hasMoreTokens()) 
vRox.setElementAt((lnteger.valueOf(t.nextToken() )), i ); 

for (s=0; s<16; s++){ 

Integer aa = (Integer)vFam.elementAt(s); 
Integer bb = (Integer)vTet.elementAt(s); 
Integer cc = (Integer)vTam.elementAt(s); 
Integer dd = (Integer)vRox.elementAt(s); 

a.addCycle(s, aa.shortValue(), bb.shortValue(), 
cc.shortValueQ, dd.shortValue() ); 

// cy = a.site[s].cycle -1 ; 

//Debug.log("Main: Site " +s+ " Cycle " +cy+ " " + 
a.site[s].dye[0].rOptic[cy]+ 

// " "+a.site[s].dye[1].rOptic[cy]+ 
// " "+a.site[s].dye[2].rOptic[cy]+ 
// " "+a.site[s].dye[3].rOptic[cy] ); 
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} 

} 

} 

catch(IOException e) { 
Debug.log("IOException"); 

} 

// Set up Melt Inverse of FAM 
for (s=0; s<16; s++){ 
for (short sec=0; sec<a.site[s].cycle; sec++) { 
//Debug.log ("Adding data to Melt " + sec + "" + 
a.site[s].dye[1 ] . rOptic[sec]); 

a.addMelt(s, sec, a.OPTICS, a.site[s].dye[1].rOptic[sec]); 
a.addMelt(s, sec, a.TEMP, (short)(60+sec)); 

} 



/* 

// Set UP for quantation. 
//100 

a.setSiteType(0, SITE_STANDARD); 
a.setConc(0, 0, 100f); 

a.setSiteType(1, SITE_STANDARD); 
a.setConc(1, 0, 100f); 

//1000 

a.setSiteType(3, SITE_STANDARD); 
a.setConc(3, 0, 1000f); 

a.setSiteType(8, SITE_STANDARD); 
a.setConc(8, 0, 1000f); 

//10 

a.setSiteType(14, SITE_STANDARD); 
a.setConc(14, 0, 10f); 

a.setSiteType(1 5, SITE_STANDARD); 
a.setConc(15, 0, 1 0f); 

// Unknowns 

a.setSiteType(2, SITE_UNKNOWN); 
a.setSiteType(4, SITE_UNKNOWN); 
a.setSiteType(5, SITEJJNKNOWN); 
a.setSiteType(6, SITE_UNKNOWN); 
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sites) 



a.setSiteType(7, SITEJJNKNOWN); 

for (int i=9; i<14; i++) 

a.setSiteType(i, SITE_UNKNOWN); 

*/ 
/* 

// Force QIC Cycle for testing 
for (int i=0; i<16; i++){ 

a.setTCycle(i, 1, (float)(10+.1*i)); 

//a.setTCyclefl, 1, 10f ); 

a.site[i].dye[1].tValid = true; 

} 

for(int i=0; io.numSites; i++) 
a.updateQuantitative(i); 

*/ 

// (site, dye, data) 

//a.dl_og(7, 1,1); //outputs threshold limits + Cycle num 

//a.dLog(7, 0, 0); // outputs data 

//a.dl_og(7, 1 , 2); // outputs raw + 2d 

//a.dl_og(7, 0, 3); // outputs threshold limits + Cycle num 

//a.dl_og(7, 0, 4); // outputs threshold limits + Cycle num + QIResult 

//a.dLog(0, 0, 5); // outputs Tlimits + TCycle num + cone (dye, all 



//a.dl_og(0, 0, 6); // outputs qtArr for given dye 
//a.dl_og(7, 1 , 7); // outputs threshold limits + Cycle num + QIC 
Cycle numbers 

//a.dl_og(7, 1 , 8); // Outputs melt data for given site. 
//a.dl_og(7, 1 , 9); // Outputs melt data peaks for given site. 

Debug \QQ^'****************^**^ M ** M *^*^*^*^**^*********** n y 
Debug. log("data4.csv, primary w Man Thresh, 
setNoiseSubtraction(true)"); 

Debug. log("setBoxCarAvg(true, 3) Quantitative "); 
Debug. log ("***************************************************^ 
a.dl_og(3, 0, 2); 

} 



llllllllllllllllllllllllllllllllllllllllllllilllllllllllllllllllllllllllll 

llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

II Used for unit testing 

void dl_og(int st, int dy, int data) { 



AttyRef#22660-0025B US 



163 



ii ii ^ 

ii ii 



int i, s, d, c; 

switch (data) { 
case 0: 
5 // data 

Debug.log("dLog: pOptic 7,* - Cy 0-44"); 

for (i=0; i<site[st].cycle; i++) 

Debug.log(" " + site[st].dye[0].pOptic[i] + 
10 " " + site[st].dye[1].pOptic[i] + 

+ site[st].dye[2].pOptic[i] + 
+ site[st].dye[3].pOptic[i] ); 
break; 
case 1 : 

^ 15 // thresh Limits, Cycle Numbers 

t! for (s=0; s<numSites; s++) 

S5 for (d=0; d<MAX_DYES; d++) 

Z Debug.logfSite " + s + 

i " Dye " + d + 

a 20 " Thresh " + getTLimit(s, d) + 

i " Cycle " + getTCycle(s, d) ); 

« break; 

fy // Prints raw + 2d data for st, dy 

*J 25 case 2: 

for (c=0; c<site[st]. cycle; C++) 
Debug.log("Site " + st + 
" Dye " + dy + 
" Cycle " + c + 

30 " raw data " + site[st].dye[dy].rOptic[c] + 

" data " + site[st].dye[dy].pOptic[c] + 
" 2D " + site[st].dye[dy].d2pOptic[c] ); 
break; 

35 // Prints dy channel TCycles and TLimits 

case 3: 
for (s=0; s<numSites; s++) 
Debug.log("Site " + s + 
" Dye " + dy + 

40 " Thresh Cycle " + getTCycle(s, dy) + 

" Thresh Limit " + getTLimit(s, dy) 

); 

break; 

45 // Prints dy channel TCycles and TLimits and Ql Results 

case 4: 



Atty Ref# 22660-0025B US 



164 



for (s=0; s<numSites; s++) 
Debug.log("Site " + s + 
" Dye " + dy + 

" Thresh Cycle " + getTCycle(s, dy) + 
5 " Thresh Limit " + getTLimit(s, dy) + 

" Result " + getQLResult(s, dy) 

); 

break; 

10 // Prints dy channel TCycles and Qn Results 

// for dye at all sites 
case 5: 
for (s=0; s<numSites; s++) 
if (useQIC) { 
15 Debug.logfSite " + s + 

" Dye " + dy + 

y " QIC Thresh Cycle " + getQICTCycle(s, dy) + 

! L i " Result " + getQTResult(s, dy) 

J ); 

O 20 } 

£ else { 

3 Debug. log("Site " + s + 

H " Dye " + dy + 

FU " Thresh Cycle " + getTCycle(s, dy) + 

l=* 25 " Result " + getQTResult(s, dy) 

S ); 
° } 

r ~" break; 

30 case 6: 

for (c=0; c<qtArr[0].length; C++) 

Debug.logf qtArr[0] Len "+ qtArr[0].length +" cone "+ 
qtArr[0][c].conc+ " Avg cy "+ qtArr[0][c].avgTCycle); 
break; 

35 

// Prints dy channel TCycles and TLimits + QIC 
case 7: 
for (s=0; s<numSites; s++) { 
for (dy=0; dy<4; dy++) { 
40 Debug.log("Site " + s + 

" Dye " + dy + 

" Thresh Cycle " + getTCycle(s, dy) + 
" QIC Thresh Cycle " + getQICTCycle(s, dy) + 
" Thresh Limit " + getTLimit(s, dy) 
45 ); 

} 
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} 

break; 

// Prints melt for given site 
case 8: 

for (c=0; c<site[st]. cycle; C++) { 
Debug.logfSite " + st + 
" sec " + c + 

" mOptic " + site[st].mOptic.get(c) + 
" d1 mOptic " + site[st].d1 mOptic.get(c) 
" Temp " + site[st].mTemp.get(c) 

); 

} 

break; 

// Prints melt Peaks for given site 
case 9: 

for (c=0; c<site[st].getMeltPeakCount(); C++) { 
Debug.log("Site " + st + 
" MeltPoint " + c + 

" dlpeak " + site[st].mPeaks[c].d1Peak 
" temp " + getMeltTemp(st, c) 

); 

} 

break; 

} 

} 

} 
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// l/lll/llll/lf/lll/l/llf/m 
II Least Squares Fit. Takes an array of points (x,y pairs) and calulates 
// the slope and offset using the 'Least Squares Fit 1 method. 

// //////////////////////////////////////////////////////////////////////// 
class LeastSquares { 

double sumX = 0.; 
double sumY = 0.; 
double sumXY = 0.; 
double sumOfXSq = 0.; 
double sumXSquared = 0.; 
int arrayLen = 0; 
double slope = 0.; 

LeastSquaresQ {}; 

// Used for quantation. 

LeastSquares(Analysis.StdElement aQ, int d) { 

arrayLen = a. length; 

for(int i = 0; i < arrayLen; i++) { 
sumX += a[i].avgTCycle; 
sumY += a[i].conc; 
sumXY += a[i].avgTCycle * a[i].conc; 
sumOfXSq += a[i].avgTCycle * a[i].avgTCycle; 

}; 

sumXSquared = sumX * sumX; 

}; 

// Used for removing background noise 
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LeastSquares(float opticQ, int start, int end) { 

arrayLen = end- start + 1; 

5 for(int i = start; i < end+1 ; i++) { 

sumX += i; 
sumY += opticfi]; 
sumXY += i * optic[i]; 
sumOfXSq += i * i; 

10 } 

sumXSquared = sumX * sumX; 

}; 

double getSlope() { 
15 if(Math.abs(sumOfXSq - sumXSquared / arrayLen) > 10E 
slope = (sumXY - (sumY * sumX / arrayLen)) / 

(sumOfXSq - (sumXSquared / arrayLen)); 

} 

else { 

20 slope = 0.; 

} 

return slope; 

} 

25 double getOffset(){ 

return (sumY / arrayLen) - (slope * sumX / arrayLen); 

} 

} 
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II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllllll 

II This object takes 2 points (x,y) pairs and calculates the slope and 
// offset. It returns the unknown (either x or y) using the equation 
// y = mx + b. 

////////////////////////////////////////////////////////////////////////// 
class LinearFit { 

double m; 
double b; 

LinearFit() {}; 

LinearFit(int x1 , double y1 , int x2, double y2) { 
m = 0.; 
b = 0.; 

if((x1 - x2) != 0) { 
m = (y1 - y2) / (x1 - x2); 
b = y1 - m * x1 ; 

} 

} 

LinearFit(float x1 , double y1 , float x2, double y2) { 
m = 0.; 
b = 0.; 

if((x1 - x2) != 0) { 
m = (y1 - y2) / (x1 - x2); 
b = y1 - m *x1; 

} 

} 
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float fitX(float x) { 
return (float) (m * x + b); 

} 

float fitY(float y) { 
if(m != 0) { 

return (float) ((y - b) / m); 

} 

else { 
return 0; 

} 

} 
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II llllllllllllllllllllllllllllllllllllllllllllllllllllllllllltllllllllllll 

II Determines the Peak and Cycle for the second derivative. It takes 3 
// points (x,y pairs) and fits a line of the 2nd order through all three 
// points, peak(y) is optic and cycle(x) is the PCR Cycle number. 

II IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIHII 
class PeakFinder { 

float peak; 
float cycle; 

double dO, d1,d2, d3; 
double r1, r2, r3; 

PeakFinder () {}; 

PeakFinder(float x1 , float y1 , float x2, float y2, float x3, float y3) { 
dO = det((x1 * x1 ), x1 , 1 , (x2 * x2), x2, 1 , (x3 * x3), x3, 1 ); 

d1 = det(y1 , x1 , 1 , y2, x2, 1 , y3, x3, 1 ); 

d2 = det((x1 * x1 ), y1 , 1 , (x2 * x2), y2, 1 , (x3 * x3), y3, 1 ); 

d3 = det((x1 * x1 ), x1 , y1 , (x2 * x2), x2, y2, (x3 * x3), x3, y3); 

if(dO != Of) { 
r1 =d1 / dO; 
r2 = d2 / dO; 
r3 = d3 / dO; 

cycle = (float) ((-1 * r2) / (2 * r1 )); 
peak = (float) (r3 - (r2 * r2) / (4 * r1 )); 

} 

else { 
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cycle = Of; 
peak = Of; 

} 

} 

\\nuiu\mnm\uuu\ninimui\niimniuiuini\uiu\iu\i 

double det(float a1 1 , float a1 2, float a1 3, float a21 , float a22, float a23, 
float a31 , float a32, float a33) { 

return ( (a11 * a22 * a33) + (a12 * a23 * a31) + (a13 * a21 * a32) - 
(a31 *a22*a13)-(a32*a23*a11)-(a33*a21 *a12)); 

} 
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